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(N ! ABSTRACT 

■ Several BL Lac objects are confirmed sources of variable and strongly Doppler- 

' boosted TeV emission produced in the nuclear portions of their relativistic jets. It 

is more than probable, that also many of the FR I radio galaxies, believed to be the 

parent population of BL Lacs, are TeV sources, for which Doppler-hidden nuclear 7- 

q ■ ray radiation may be only too weak to be directly observed. Here we show, however, 

5— 1 ' 

' that about one percent of the total time-averaged TeV radiation produced by the ac- 

tive nuclei of low-power FR I radio sources is inevitably absorbed and re-processed 
by photon-photon annihilation on the starlight photon field, and the following emis- 
s_j ■ sion of the created and quickly isotropized electron-positron pairs. In the case of 

$ ■ 

the radio galaxy Centaurus A, we found that the discussed mechanism can give a 
distinctive observable feature in the form of an isotropic 7-ray halo. It results from 
the electron-positron pairs injected to the interstellar medium of the inner parts of 
the elliptical host by the absorption process, and upscattering starlight radiation via 
the inverse-Compton process mainly to the GeV— TeV photon energy range. Such a 
galactic 7-ray halo is expected to possess a characteristic spectrum peaking at ~ 0.1 
TeV photon energies, and the photon flux strong enough to be detected by modern 
Cherenkov Telescopes and, in the future, by GLAST. These findings should apply as 
well to the other nearby FR I sources. 

Key words: radiation mechanisms: non-thermal — gamma-rays: theory — galaxies: 
active — galaxies: jets — galaxies: individual (Centaurus A) 
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1 INTRODUCTION 



Cent aurus A is the most nearby active galaxy 1 , hosted by a powerful elliptical (see a monograph 



by 



Israel 



199a) . At low frequencies it reveals a giant (8 deg x4 deg or 500 kpc x 250 kpc) and 
complex radio structure, with a total 5 GHz energy flux of 3.4 x 10~ n erg cm~ 2 s^ 1 , containing 
a one-sided kpc-scale jet, few-kpc-long inner lobes, and extended outer lobes. The host galaxy 
is of 14' x 18' optical size, has a blue luminosity of 7.5 x 10 10 L Q , and contains a ~ 10 8 M e 



supermassive black hole in the center (IMarconi et al 



2006). The famous 'dark lane' pronounced 



within the galactic body, being in fact an edge-on disk of rotating gas, dust and young stars, is 
most probably a remnant of the merger with a spiral galaxy, which happened some 10 s — 10 9 years 



ago. The kpc-scal e jet and the active nuc 



recent analysis by 



Hardcast le et al J 2003 ; 



eus are confirmed s ources of X-ray emission (see most 

II 1 

Kataok a et al.1120 06. and references therein). 



Cen A source is classified as an FR I radio galaxy, although some of its morphological char- 
acteristics (double-double structure, one-sided jet) are not typica l for this class of AGN s. FR I 
galaxies are believed to be a parent population of BL Lac objects ( Urrv & Padovanil 19951) . some 
of which are confirmed sources of very high energy (VHE) 7-ray emission. In fact, the nuclear 
part of the Cen A jet was modeled in terms of a 'misaligned' (and therefore Doppler-hidden) BL 



Lac, since the jet viewing angl e in this source is expected to be much 



clination of the blazar jet (e.g., 



Steinle et al. 1998 



Baile y et al. 



1986; 



Mor ganti et al. 



1991; 



arger than the typical in 



Botti & Abraham 



1993 



Chiaber ge et al.H2001h . As such, Cen A was considered for some time as a po- 



tential source of TeV radiation, for which jet misalignment effects are compensated by the source's 
proximity. This expectation was strengthened by the positive detection of the considered object at 
Me V— GeV energy range by all the instruments on board CGRO during the period 1991-1995 



(Steinle et al 



1998, and references therein). Nevertheless, the observations perf ormed till now at 



2005) 



the TeV energy range result ed in the upper limit for such an emission, only (Ro well et al 
Aharonian etal]l2004 



1999; 



In this paper we investigate the emission resulting from re-procession of the potential 7-ray 
flux of the active nucleus in the Cen A radio galaxy via annihilation of the nuclear 7-rays on the 
starlight photon field due to the host galaxy, followed by the synchrotron and inverse-Compton 
cooling of the created electron-positron pairs. We show, that the analysis of this mechanism — 
which also applies to other FR I sources — can give us interesting constraints on the unknown pa- 
rameters of the Doppler-hidden and heavily obscured Cen A nucleus, as well as of its host galaxy 



1 Distance d = 3.4 Mpc; 1 arcsec corresponds to 16 pc. 
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and outer parts of its radio outflow. We note that 7-ray opacity in blazar sources was considered be- 
fore only in the co ntext of nuclear target radiation fields, due to accretion discs, broad-line regions . 



or dusty tori (e.g., 



Wag ner et al. 



1995; 



Sikora et al 



1994 



Bottcher et al 



1997; Wang 



1994; 



Rlandford & Levinson 



1995 



Dernier & Schlickeiser 

200oJ). In the case of the FR I/BL Lac sources, 



however, such radiation fields are relatively weak, or even absent. Thus, the main source of the 
opacity for the TeV photons emitted from their active nuclei is the starlight radiation 2 . This radia- 



are ho sted by giant elliptical galaxies (e.g., 



Colina & de Juan 


1995; 


Urrvetal. 


2000; 


Heidtetal. 



2004), for which spectral and spatial distribut ion of the ste 



(see in this context discussion in 



Stawarz et al. 2005 



2006a||bJ). 



lar emission is relatively well known 



Below, in section 2, we present details regarding the calculation of the optical depth for photon- 
photon pair production by the Cen A nuclear 7-ray emission on the starlight photon field of the 
elliptical galaxy, and the energetics of the created electron-positron pairs. In section 3, we discuss 
further evolution of such particles injected into the interstellar medium of the elliptical host, cal- 
culating in particular the resulting synchrotron and inverse-Compton fluxes. Section 4 contains a 
final discussion and conclusions. 



2 CALCULATIONS 



2.1 Opacity 



In order to calculate the opacity for the 7-ray photon beam propagating through the galaxy, one 
has to specify the spectral and spatial distribution of the stellar photon field(s). Starlight surface 



brightness in elliptica 



galaxies, includin g those with active nuclei, are typicaly well fitted by an 



empirical Nuker law (ILauer et al.l ll995). In terms of monochromatic starlight intensity, this law 
can be expressed as 

a-\ —(b—d)/a 

(1) 



I(r) = I h 2 {b - d)/a 



-d r 



1 + 



where r is the distance from the galactic nucleus, is the 'break' radius, and lb = /(rb). Such 
profiles imply a power-law dependance J(r) oc r~ d for r < rb, and I(r) oc r~ b at larger distances. 



We note, that ellipticals h osting radio sources may contain relatively large amounts of cold dust, prounounced at far infrared frequencies 

for which the exact location is unknown. As mentioned before, the Cen A system itself is bright in 



Golombek t 



et alj|l988l Knapp et alj|l99pl) 



infrared frequencies due to the dust and young stars emission, restricted however to the 'dust lane', i.e. relatively thin but extended disc-like fea- 
ture (roughly perpendicular to the jet axis), being a remnant of the spiral which merged recently with the elliptical host. A role of this feature in 
absorbing and reprocessing the nuclear 7-rays is briefly discussed in Appendix A. 



Stawarz et al. 



de Ruiter et al. 



(2005) found that in the case of weak radio galaxies (selected from the B2 sample 
with the criterium of not showing significant dust emission), the host galaxies are characterized by 
the typical values of a = 1.9 , b = 1.6, and d = 0.02. In a specific case of the Cen A host galaxy, 
Capetti & Balmaverdel J2005I) give a = 1.68,6 = 1.3, d = 0.1, and a break radius r b = 2.56" = 41 
pc. These values are considered hereafter, allowing to approximate the starlight emissivity, j(r) oc 
r" 1 J(r), in the NGC 5128 (Cen A host) galaxy, as 



(2) 



where £ = r/rv Here e = e/m c c 2 is the starlight photon energy e in m c c 2 units, jy is the total V- 
band galactic emissivity, g(e) describes the V-band normalized spectral distribution of the stellar 
photon field, and 

1 Kr) 



-0.7143 



M.0 = - I ^ = 1.64 r 1 '' (1 + £"*) (3) 

r I — II — I 

is the radial depedance function. The considered galaxy is 14' x 18' in optical size (Israel 1998), 
and hence we take the average terminal galactic radius of 16', or £ t = 7"t/?"b = 375, leading to the 
galactic volume (V) integral 

fit 



n 



h(£) dV = 4nrl I ( z h{i) d£ = 1.83 x 10 3 r 



3 ^3 

b 



(4) 



'V J0 

(assuming spherical symmetry). For the spectral distribution of the N GC 5128 starlight, we assume 



Silva et al 



(1998), restricted 



the template spectrum of a powerful elliptical galaxy as provided by 
to the photon energy range between e min = 10~ 7 and e max = 10 -5 . This restriction implies that we 
consider only the direct stellar emission of the evolved red giants (constituting the main body of 
the elliptical host) and their winds, but not the far infrared emission resulting from the reprocession 
of the starlight photons by the cold galactic dust (see Appendix A). We normalize it to the V-band 
radiation, so that g(e = h/m c c\y) = 1, where Ay = 0.55 u m. The resu l ting s pectral function g(e) 



Silva et al. 



is shown in Figure 1 (top-panel): values read directly from 
crosses, and the broken power-law approximation considered hereafter, 

' 2.67 x 10 14 e 2 - 44 
9.25 x 10^ 4 e~ a57 



( 1998) are represented by 



9\t) 



< 



for 10" 7 - ^ e ^ 10~ 5 - 8 
for 10" 5 - 8 < e < 10" 5 - 3 



(5) 



3.54 x 10~ 31 e" 5 74 for 10" 5 - 3 < e < 10" 5 - 



is shown as a solid line. 



The apparent V"-band magnitude of NGC 5128 is m v = 6.98 (Isr aeli 1 199 81) . This gives the 
monochromatic V-band galactic luminosity 
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Figure 1. Top: Spectral energy distribution of the starligh emission of a template giant elliptical; values read directly from Sil va et all 1 1998t) 
are represented by crosses, and the broken power-law approximation introduced in this paper is shown as a solid line. Bottom: Optical depth for 
annihilation of the nuclear 7-ray emission on the starlight photon field. 



log(^V) = 50 - 078 + 2 log (ttM -°- 4 K- M)~ c v = 43.82 , (6) 
\erg/s/ \Mpc/ 

where the distance of Cen A is = 3.4 Mpc, the extinction is Ay = 0.381, and cy = 4.68. Since, 

by the definition, 

L v = 4vr [ [ej e (0] Y dV = 47rj v H , (7) 
Jv 

one obtains jy = 1.42 x 10~ 21 erg cm -3 ster -1 s -1 . Note, that by means of standard radiative 
transfer formulae, the spectral intensity within a solid angle fi, which can be obtained by inte- 
grating spectral emissivity along a ray as = J dl, is related to the differential photon 
number density, n e (£, Vt), through the expression I u (Cl) = chen e (^, Q). This corresponds in fact 
to the neglegible absorption of the starlight along the light path, the approximation which may 
be considered as being in conflict with the presence of the dust lane in Cyg A host. It is however 



6 Stawarz et al. 

good enough for the purpose of the presented evaluation, since integration over the whole extended 
elliptical volume reduces the starlight absorption effects. As a result, 



e- 2 



n £ (£,0) = / [ej e (Z)]dl . (8) 



m e c 3 

The optical depth for photon-photon annihilation, computed for the case of a monodirectional 
beam of 7-ray photons with a dimensionless energy e 7 , propagating through the stellar photon 
field of the host galaxy from the active center up to the terminal distance r t , is 

r(e 7 ) = / dr dn(l — zu) cx 77 , (9) 



where zu is the cos function of the angle between the 7-ray photon and the incident starlight 
photon, dn = n e (£, f2) de dVt is the differential starlight photon number density, and 



cr 77 (e 7 , e, zu) = (1 - (3 2 ) 
16 



(3-/3 4 ) ln(i±|) -2/3(2-/3 



is the photon-photon annihilation cross section (iGould & Schredei 



(10) 



1967), where 



/ 2 \ 1/2 

/3= 1 n r , (ID 

is the velocity of the created electron/positron (computed via the conservation of four-momentum) 
in the appropriate center-of-momentum frame. By choosing dVL = dcf) dvu, the differential starlight 
photon number density (see equation 8) reads as 

£ — 2 ~, /''/max 

n t ^,Q) = £ / [eje(C)]^ , (12) 

m eC Jo 

where rj = l/r^, ( = a/£ 2 + r\ 2 — 2^rjzu, and the integration upper limit is r] max = £cu + 
\f£, 2 zu 2 — i 2 + £ t 2 . This gives finaly 

r(e 7 ) = ^ / dzu(l—zu) / rfe<T 77 (e 7 , e, tu) — — / d£ / d?7 /i(C) , (13) 

where the threshold energy is e t hrc = 2/e 7 (1 — w). This optical depth as a function of the 7-ray 
photon energy is shown in Figure 1 (bottom-panel). As can be seen, it possesses a broad maximum 
at e 7 w 1 - 10 TeV, with the maximum value r(2 TeV) w 0.0145. 

Let us fix the 7-ray photon energy at e 7 = 10 6 6 , and compute again the appropriate optical 
depth for photon-photon annihilation as a function of the upper bound in integration over the 
distance from the core, r = r(£), replacing simply constant £ t in equation 13 with variable £. 
Such a profile is shown in Figure 2 (solid line): as can be seen, only about 40% (r w 0.006) of 
the annihilated radiation is in fact absorbed within the innermost galactic volume of radius r b , and 
the remaining 60% is absorbed at distances 77, < r < 100 rb. The obtained profile for r can be 
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Figure 2. Left: Optical depth for 2 TeV photons as a function of the distance from the active nucleus (solid line). Right: Profile of the energy denisty 
of the starlight photon field (dotted line). 

compared with the spatial distribution of the galactic radiation. Noting the definition for the photon 
field spectral energy density, 

U u = — I v (n)dm , (14) 

c J-l 

one may find the appropriate bolometric energy density profile for the starlight emission as 

2irr h jy f +1 f mM 

U rad (0 = /boi [vU v ] y = /boi dzu drj h(() , (15) 

c J -1 Jo 

where / bo i = 2.5 is the V-band bolometric correction for the adopted here stellar spectrum (equa- 
tion 5). This profile is shown on Figure 2 as a dotted line. 

The spectral and spatial shapes of the optical depth for the photon-photon annihilation dis- 
cussed above for the Cen A radio galaxy is universal, and should hold also for the other low-power 
radio sources hosted by elliptical galaxies. That is caused by the appropriate scaling r oc jy 
(equation 13), jy oc Ly/H (equation 7), and H = nr\ (equation 4), where k = K,(a,b,d) is a 
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numerical factor depending on the parameters of the Nuker profile (for example, k = 1.83 x 10 3 in 
the case of NGC 5128, as given in the equation 3). These relations gi ve r oc Ly / 'nr h . Since the pa- 
rameters of the NGC 5128 Nuker profile (a = 1.68, b = 1.3, d = 0.1 Icanetti & Ba1maverdJ20 05^ 
are s imilar to the average hosts of B2 radio sources (a = 1.9, b = 1.6, d = 0.02 Ide Ruiter et all 



2005), since the spectral energy distribution of the starlight em ission considered here is based on 



a universal spectrum of the elliptical galaxy (iSilva et al 
linearily with the V-band galactic luminosity, 



1998), and since the break radius scales 



_Tb_ 



A' 



( de Ruiter et al 



15 



org/ s 



(16) 



2005), the optical depth for the 7-ray photons shown on Figures 1 and 2 for the 



case of Cen A radio galaxy indeed should apply also to the other analogous sources. We note for 
example, that with the value of L v found previously for the NGC 5128 galaxy, the simple scaling 
introduced in equati on 16 implies the break radiu s ~ 66 pc, which is very close to the observed 
value of r b « 41 pc ( Capetti & BalmaverdJ2005 ). On the other hand, the starlight energy density 
profile [/ ra d oc jy rb oc r^ 1 is supposed to change for different sizes (and hence luminosities) of 
the elliptical host when compared with the particular profile shown on Figure 2 for the case of 
NGC 5128 parameters. 



2.2 Energetics 

The obvious source of 7-ray photons in the Cen A system is the active nucleus (containing a 
relativistic jet), which is most probably responsible for production of the flux detected by OSSE 
(0.05 - 4 MeV), COMPTEL (0.75 - 30 MeV) and EGRET (0.1 - 1.0 GeV) instruments on board 
CGRO i n the period 1991 - 1995, with the total observed 50 keV - 1 GeV luminosity ~ 3 x 10 42 



erg s _1 (Ste inle et al.l Il998. and references therein). This emission peaks at « 0.1 MeV, with 
the maximum energy flux about ~ 10~ 9 erg cm~ 2 s _1 . At a very high 7-ray energy range, a 3a 
detection of C en A with the p hoton flux F(> 0.3 TeV) = 4.4 x 10~ n ph cm -2 s _1 was reported in 



the seventies ( Grind lay et al 



1975|) . Subsequent CANGAROO observations resulted in the upper 



limit F(> 1.5 TeV) < 5.45 x 10 12 ph cm 2 s 1 for the point source centered at the Cen A 
nucleus, and F(> 1 .5 TeV) < 1.28 x 10 -11 ph cm" 2 s _1 for the extended region with radius 14' 
( Rowellet all 1999 ). The most recent HES S observations gave F( > 0.19 TeV) < 5.68 x 10~ 12 
cm -2 s -1 for the point-like active nucleus ( Aharonian et al.ll2005 ). 



ph 



Here we are interested in the total, time-averaged and angle- averaged (i.e., 'calorimetric') flux 
of TeV photons produced by the active nucleus and 'injected' into the host galaxy as well as into 
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the large-scale radio structure. With the most recent HESS results, we can only put some upper 
limits for it. Namely, assuming a standard apectral index a 7 = 1 (equivalent to the photon index 



+ 1 



2) at the considered photon energy range, the relation between the integrated 



photon flux and the monochromatic flux energy density at any photon energy e ^ e is simply 
[£S e ] = [e F(> e )]. This flux is related to the emitting fluid (jet) intrinsic monochromatic power 
(assumed to be isotropic in the jet comoving frame) radiated in a given direction, dV / dQ! = 
L'/An, by the expression 



Ws v ] 



80! 



1 *Lc # 

£ T nur 4tt 



(17) 



where r nuc and 5 nnc = T n ^ c ^1 — a/1 — r nu 2 c cos 9 \ are, respectively, L orentz and Doppler 



factors of the nuclear portion of the jet, and 9 is the jet viewing angle (e.g., 



Stawarz et. al 



Sikora et. al 



1997; 



2003). On the other hand, the total power radiated into the ambient medium being 



of interest here, is 



A 3 Fir' 1 

■Triiir. 2 



5zt sin 9 d9 = U 



and hence 



47r4r nuc 5 nu 3 c [e F(> e Q )\ 



(18) 



(19) 



With the HESS photon flux F(> 0.19 TeV) < 5.68 x 1(T 12 ph cm" 2 s" 1 (corresponding to the 
few-arcmin-integration area centered on the Cen A nucleus) one obtains the upper limit for the 
total injected monochromatic power L in j < 2.4 x 10 39 r nuc 5~ u 3 c er g s" 1 . With the prefered values 



9 ~ 50 deg —80 deg inferred from the VLBI radio observations dJones et al 



1996; 



Tingav et al 



19981) . and r nuc ~ 10 widely considered as a typical value for the bulk Lorentz factor of sub- 
parsec scale AGN jets, this reads as L in j < 10 42 — 10 43 erg s -1 . Below we take conservatively 
L in j = 10 42 erg s _1 as an upper limit, noting that such a luminosity would correspond to less than 
10% of the minimum total kinetic power of the Cen A jet, estimated to be L j > 10 43 erg s _1 



Kraft^Ui 



2003). As shown in the 



from the radio/X-ray lobes' energetic (see Clarke 
previous section, about 1% of this power is absorbed on the starlight photon field within the central 
100 r b ~ 4 kpc region of the host galaxy, and is thus converted to an electron-positron population, 
mainly in the 0.1 — 1 TeV energy range. 

The nuclear 7-ray emission postulated here is expected to be Doppler-boosted within the nar- 
row cone characterized by the opening angle T~^ c < 6 deg. Therefore, it is strongly Doppler hidden 
when viewed from 9 ^ 50 deg (see equation 17). If the observer was located within the beaming 
cone of this emission however, he would detect a flux corresponding to the isotropic luminosity 
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Figure 3. Different approximations for the pair injection energy spectrum (corresponding to the photon flux of the primary -y-ray s oc e^ 2 ) in nor- 
malized units: (i) Q(y) oc y~ 2 t(2j) (open circles), (ii) Q(7) calculated using the pair production rate expression introduced by Aharonian et al. 
1 1983, solid line), and (iii) broken power-law approximation introduced in this paper (dotted line). 



L(0) = (5l uc 0=o /r nuc ) L' ps r^ uc V < 10 44 erg s \ Such values are consistent w ith luminosities 



Kata rzyhski et al.11200 6). whi ch are be 



observed from the TeV-detected BL Lac objects (see, e.g., 
Keyed to be beamed analogues of the low-power FR I radio galaxies like Cen A ( Urn & Padovani 



1995). 



In order to find the energy spectrum of the created electrons and positrons, we note that the 
pair production rate for small values of the optical depth can be estimated as 

dr(e 7 )\ 

J. ,^-= 2 ^ 



Q(7,r) oc <^ n 7 (e 7 ) x 



dr 



(e.g. 



Coppi & Blandlor' 



(20) 



1990; 



Bottcher & Schlickeiser 



1997), where 7 is the Lorentz factor of 



the created particles, and n 7 (e 7 ) is the photon spectrum of the primary 7-ray photons. For the 
latter we assume power-law form ra 7 (e 7 ) oc e 7 7 , where T 7 = a 7 + 1 is the photon index. When 
averaged over the galactic radius, the pair production rate therefore scales as 



Q( 7 )oc 7^(27) 



(21) 
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This function is shown in Figure 3 in normalized units (open circles) for T 7 = 2, taken in this paper 
as a typical photon index of blazars' TeV emission. As illustrated, Q{^) is peaked at 7 ~ 10 5 5 , 
and decreases as oc 7~( r T +0 - 5 ) for 7 > 10 6 . Such behavior is in fact expected, since for e 7 > 10 6 6 
the optical depth evaluated previously is roughly r(e 7 ) oc e 7 0,5 . Figur e 3 shows also the pair 



injection energy distribution calculated using the expression introduced by 



Aharonianet al. (1983 



solid line), integrated over the spatial-averaged spectrum of the starlight photons (equation 12; see 
also equation 34 bel ow) and the assumed ph oton spectrum of primary 7-rays, oc e~ 2 . We note that 
the approximation of lAharonian et al . ( 1983) for the energy distrib ution of the secondary electron s 
provides accuracy better than few% (see detailed calculations in 



Bottcher & Schlickeiser 



1997). 



Sma ll differences between the two presented estimates come from the fact that the function given 



by 



Aharonian et al. 



( 19831) corresponds to the isotropic distribution of the soft photons, while 
the optical depth calculated by us (i.e. equations 20-21) includes the anisotropy of the starlight 
radiation field. Finally, Figure 3 shows also the simplest broken power-law approximation for 
the injection pair spectrum considered hereafter for the purpose of the following calculations: 

Q(7) oc const for 70 = 10 5 ^ 7 ^ 7b r = 10 6 , and Q{^) oc 7 -2 - 5 for 7 > 7b r - We normalize it to 



the monochromatic 1 TeV power L in j specified above, to obtain 



Q{i)= I^fJ where Z( 7 ) 

7br m eC 2 V ga l 

and zero otherwise. Here r = r(e 7 = 
pair injection is taking place. 



(22) 



1 for 70 ^ 7 ^ 7br 

(7br/7) 2 ' 5 for 7 > 7b r 
1 TeV) Rj 0.01, and V ga i is the galactic volume to which the 



3 IMPLICATIONS 

The electrons and positrons produced by annihilation of the VHE 7-ray nuclear photons on starlight 
radiation are expected to have an initial spectrum peaked at 0.1 — 1 TeV energies. After being in- 
jected into the galactic medium, they are quickly isotropized by the ambient magnetic field, and 
radiate via the synchrotron and inverse-Compton processes. The evolution of these electrons and 
details of their emission spectra depend therefore on the properties of the interstellar medium of 
the elliptical host. 3 

3 We also note, that since the nuclear 7-ray emission is expected to be beamed within the narrow cone with an opening angle of less then 10 deg, as 
discussed in the previous section, some part of it illuminates also the large-scale radio outflow, i.e. the kpc-scale jet (extending in the case of Cen A 
from ~ 0.2 kpc up to ~ 4 kpc from the center), and the 'inner' lobe (~ 4 — 5 kpc; see lfiurns et alil983llciarke et alJl99 j) . The evolution of 
the electron-positron pairs injected via the photon-photon annihilation into these regions can differ substantially from the evolution of the particles 
injected into the body of the elliptical host as considered below. This problem will be investigated in a subsequent paper. 
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3.1 Elliptical host 

Moss & Shukurovl (|1996|) argued, that elliptical galaxies have no ordered large-scale magnetic 
field, but only an unresolved random component. They further suggested, that the latter is due to a 
'fluctuation dynamo' driven by the turbulent motions of the interstellar matter caused by type I su- 
pernovae and stellar motions. The latter ones are expected to be characterized by the Kolmogorov- 
like energy spectrum and an injection scale ~ 3 pc, while the former ones by a steeper (shock-like) 
spectrum and larger injection scale, ~ 300 pc. The resulting random magnetic field is expected to 



be characterized by the average ga 
and the correlation scale ~ 100 pc 



actic intensity ~ 3 uG (r eaching ~ 10 /iG in the central parts) 



Moss & Shu kurovl (|1996) demonstrated that their model is con- 



sistent with all the observational constraints (in particular with the depolarization studies). In the 
case of Centaurus A, and in general all the ellipticals hosting radio-loud AGNs, a galactic magnetic 
field can be even higher than this, especially close to the galactic center, possessing even some reg- 
ular component due to polution of the interstellar medium by the magnetized plasma transported 
from the active core in the form of the jets. However, for the purpose of order-of-magnitude eval- 
uations, below we take the characteristic values for the NGC 5128 elliptical host's magnetic field 
-Bgai ~ 3 — 10 [iG, assuming that it consists solely of the (Alfvenic) turbulent component with 
the maximum wavelength A max ~ 100 pc and Kolmogorov energy spectrum W(k) oc k~ q , where 
q = 5/3. 

With the interstellar medium parameters as discussed above, the mean free path of the created 



elect ron-positron pairs for resonant interactions with the turbulent Alfven modes, is (Sch lickeisei 



1989) 



A P 



A, 



f 3 A^~ 0.82X7^5:^ pc 



(23) 



where 7 6 = 7/IO 6 , r g = >y m e c 2 /eB gal ~ 5.5 x 10 5 7 6 B_\ pc is the electrons' gyrora- 
dius, and = B^/10 uG (as a general reference for the particle-turbulent wave interac- 



tions see 



Schlickeiser 



2002). Such interactions lead to quick isotropisation of the injected elec- 



1/3 R -i/3 



trons. In particular, the appropriate isotropisation time scale, t iso ~ 3A /c ~ 8 7g /0 -ET 5 '~ yrs, 
is a few orders of magnitude shorter than the radiative cooling time scale (see below). The dif- 
fusive escape time scale from the central parts of the elliptical host (R ~ 100 r b ~ 4 kpc) is 
t esc ~ 3 R 2 /X e c ~ 2 x 10 8 7g^ 3 B 1 /^ yrs. We note, that the re-acceleration of the radiating par- 
ticles by resonant scattering on the Alfven modes is not expected to be efficient enough to keep 



the electrons around 7 ~ 10 6 . In particular, the time scale for this process, t acc 
4.5 x 10 6 7g^ 3 BZl^ 3 yrs, is much longer than the radiative cooling time scale. In the above, 
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Figure 4. Time scales for the radiative cooling (dark gray), turbulent re-acceleration (light gray) and diffusive escape (gray) of the electrons injected 
into the interstellar medium of the elliptical host, for the range of the galactic magnetic field 3 — 10 fj,G. The values corresponding to -B ga i = 10 
fj,G correspond to the lower bounds for the radiative losses and the acceleration time scales, and to the upper bound for the escape time scales, 
respectively. 



v A = Pac « B gal (4:7rm p n 



gas ) 



-1/2 



~ 10 3 5_5 c is the Alfven velocity expected for the av- 



erage number density of the cold gas within central (< 4 kpc) parts of the NGC 5 128 galaxy 



n 



3x10 3 cm 3 . In this context we note that the most recent analysis presented by lKraft et al 



(2003) indicates that the hot gaseous component of the galaxy discussed here is characterized by 
a tempreature kT ~ 0.3 keV and a central number density ~ 4x 10~ 2 cm -3 , which is roughly 
constant within the radius r ~ 0.5 kpc and decreases further away as oc r~ 12 . Such a behavior 



is consistent with general properties of the elliptical galaxies (Mat hews & B righent 



2003). Thus, 



one can conclude that the TeV energy electrons injected into the interstellar medium (via annihi- 
lation of the 7-ray emission of the active center) are effectively confined to the elliptical body and 
quickly isotropized by the galactic magnetic field, and therefore radiate all their energy there by 
inverse-Compton upscattering of the starlight photons and the synchrotron process (see below), 
before being re-accelerated by the turbulent processes. As a result, one should expect formation of 
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Aharonian et al. 



(1994), 



a small-scale (galactic) version of the 'isotropic pair halos' discussed by 
restricted however to the first generation of the secondary photons. 

Below we investigate in more detail the emission spectrum of electron-positron pairs created 
within the centaral 100 r b parts of the elliptical host's interstellar medium due to annihilation of the 
nuclear VHE 7-rays on the starlight photon field. First, we introduce the averaged energy density 
of the starlight photon field, assumed in this section to be isotropic within 100 r b , 

eU e « ([vU„] y )e x g(e) (24) 



(see equations 5 and 15), where 

MJma; 



(WUviv/t = r £rr 7, = 1 d£ dm drj h(Q , (25) 

Jo J-i Jo 



It** C ^ c 

and r cr = £ cr r b is some particular critical radius over which the averaging is performed. The total 
averaged starlight energy density is 

/*£max 

(Orad) = / U e de = / bol ([l/tf„] v )f . (26) 

With £ cr = 100 and other parameters as given above, one obtains {[fU v \ v )^ 10~ n erg cm -3 , 
and (Z7 ra d) ~ 2.5 x 10~ n erg cm -3 , which is still much higher than the energy density of the 
galactic magnetic field, {7 B = B^ al /8n k4x 10~ 12 B^ 5 erg cm -3 . The time-scale for the radiative 
losses of the electrons with Lorentz factor 7, including synchrotron and inverse-Compton losses 
(in both Thomson and Klein-Nishina regimes) can be then estimated as 
. . 3 m e c 

*rad(7) = "j 77 /-1 , j-i \ , (27) 

where q = (U iad )/U B ~ 7.B~J, and 

Fkn = 7fh I n , r ^ 1 ■ 5 d€ (28) 
Uwd) A (1 H- 47e) 



v - 



(Moderski et al 



2005). The simplified form of the function Fkn introduced above is in fact a very 



accurate approximation for the values 4 7 e < 4 7 in j e max < 100, as considered in this paper (with 
e max = 10~ 5 ). The time scale £ ra d(7)> together with the turbulent re-acceleration time scale and the 
escape time scale, is shown in Figue 4 as a function of the electron Lorentz factor 7 for the range 
of £? ga i = 3 — 10 fxG (note that ^0(7), not shown in this figure, is much shorter than any other of 
the time scales plotted). As can be seen, t esc 3> max(t acc , £ rad ) for all 7 < 7 br , and t acc 3> t ra d for 
7 > 70. All the electrons with Lorentz factors 7 < 7kn = l/4e max ~ 10 4 5 cool mainly via the 
inverse-Compton losses in the Thomson regime. Meanwhile, the electrons with 7kn < 7 < 7cr, 
where qF KN (j CI ) = L lose their energy mainly via the inverse-Compton emission in the Klein- 
Nishina regime. In the case of B sa \ = 3 /xG, one has 7 cr ~ 7 br , while for B ga \ = 10 jiG the 
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Figure 5. Energy distribution of the electrons injected to the interstellar medium of the elliptical host, for the galactic magnetic field 3 fiG (solid 
line) and 10 fiG (dashed line), in normalized units. Shaded region indicate the electrons' energy range in which neglected re-acceleration effects 
are expected to become important. 

critical electron Lorentz factor 7 cr is a factor of a few lower than that. Finally, the main process 
responsible for the cooling of the electrons with 7 > 7 cr is synchrotron radiation. 

The resulting electron energy distribution, n c (^), ignoring re- acceleration and escape effects, 
can be found from the continuity equation 

<9n e (7) _ d_ 
dt cfy 

where (3(7) denotes injection of high-energy electrons through photon-photon annihilation, and 



{|7|cooi^(7)} + Q(7) , (29) 



1 7 1 cool — |7|syn + \i\ic is the total rate of the radiative cooling. The standard formulae give 

4 CO"t 2 4 CO"t 2 

syn = o ^b7 , and 7 ic = 9^kn (30) 



(Moderski et al 



2005). Note, that for the parameters considered here, and electron Lorentz fac- 
tor 7 = 10 6 , the relative importance of the inverse-Compton and synchrotron energy losses, 
|7|ic/|7|syn = qF kn , is roughly ~ 3 for B gal = 3 /xG, and ~ 0.3 for 10 /iG. Thus, the energy 
injected into the created pairs is re-radiated via their synchrotron and inverse-Compton (in the 
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Klein-Nishina regime) processes in roughly comparable amounts. As for the function Q{^), we 
use approximation for the electrons freshly injected to the galactic volume V ga i introduced by 
equation 22, to obtain the steady-state solution 

S,diQ{i) 3rL inj j.dilH) 

UAl) 4(7 T T 2[/ B (l + gF KN ) 4 C <7 TT 2 r y galT 2[/ B (l+ g FKN) ' 

This spectrum, in normalized units, is shown in Figure 5 for £> ga i = 3 /iG (solid line) and 10 /iG 
(dashed line). The shaded region indicates the energy range for which the turbulent re- acceleration 
effects (neglected here) are expected to become important. As can be noted, at low electron en- 
ergies 7 < 7kn, for which the inverse-Compton cooling in the Thomson regime dominates, the 
energy spectrum has a standard form n e (7) oc 7~ p with p = 2, as expected in the case of the con- 
tinuous injection of flat-spectrum (Q(7 < 7t>r) oc const) particles, followed by the |7|i C /T oc 7 2 



energy losses. However, for 7 > 7kn the electron spectrum flattens, as a resu 
nant inverse-Compton/Klein-Nishina cooling (see the most recent discussion in 



2005 



t of the domi- 



Moderski et al 



and references therein), although in the case of B ga \ = 10 /iG this effect is relatively weak 
(the effective spectral index p ~ 1.6, to be compared with p ~ 1.4 for £> ga i = 3 /xG). Mean- 
while, at 7 > 7 cr the electron sepectral index in both cases increases to p = 3.5, as expected 
for the dominant synchrotron cooling |7| syn oc 7 2 of the continuously injected steep-spectrum 
(Q(l > 7br) oc 7~ 2 ' 5 ) electrons. 

With the evaluated electron energy distribution, one can find the energy spectrum of the result- 
ing synchrotron and inverse-Compton emissions of the created pairs. In general, for the isotropic 
distribution of seed photons and particles, the respective luminosities can be simply written as 
[eL £ ] syn/ic = 4vr V ga i [ej e ] syn/ic , where 

(32) 



1 ^(7)7 2 

[ e JeJsyn ^ . |7|syn m e c 



7=V( 3 / 4 ) «>yn (B cr /B) 



(using 5- approximation for the synchrotron emissivity; see 

2vr m 2 c c 3 /h e = 4.4 x 10 13 G, and 

3a T m c c 3 2 f f J" iso (7, e, e ic ) 

— e ic / c?7 / de n e (7)n rad (e) 

16vr J J 67 



Crusius & Schlickeiser 



1986), B cr = 



[ej e J ic = — — — e ic J d'y j de n c (7)n rad (e) — . (33) 



Here e syn / ic denotes now the dimensionless energy of the synchrotron/inverse-Compton photons, 
electron energy distribution n e (j) is given by the equation 31, seed photons' spectrum 

n rad (e) = r = YT~ 6 #( e ) (34) 

is specified by the equation 5, and finally 



, 2 . (icv)Hi-v) 



^i S o(7, e, e ic ) = 2 V \nV + V + 1 - 2V 2 + 2 ^ +JCV) (35) 
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Figure 6. Synchrotron and inverse-Compton spectra of the electrons injected into the interstellar medium of the elliptical host, for the galactic 
magnetic field 3 fj,G (solid lines) and 10 fj,G (dashed lines), and the injection luminosity L; n j = 10 42 erg s — x . Dotted lines indicates the 
0.05 — 1000 MeV spect rum of the Cen A source as constrained by different instruments on board CGRO (1991-1995; 'low' and 'intermediate' 
states ; ISteinle et all 19981) . 



with 

/C = 4 e 7 and V = 



4f 7(7 ~ ftp) 



where 1/ ! ' P £ 1 dBlumenthal & GoulJl97 ()l. This leads to 



TlW) 



2 74 1 + qFi 



KN 



and 

[eL, 



9qrL i 



mj 2 



7=v /(3/4) £syn (Bd/B) 



min(e n 



max(e min ,e low ) 6 3 7 4 (1 + ? F KN ) 



d 7 'X(7') 



(36) 



(37) 



(38) 



16/bol7br 

where ei ow = £10/47(7 — e ic ), and e up = £^7/(7 — e ic ). The evaluated luminosities are shown 
in Figure 6, for L in j = 10 42 erg s _1 and the galactic magnetic fiel £> gal = 3 /iG (solid lines) 
and 10 fiG (dashed lines). In the case of the inverse-Compton emission, spectral index in the 
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photon energy range e ic = 10 3 — 10 6 (corresponding roughly to the electrons energies 7kn < 
7 < 7brX is ttk/KN ~ 0.36 for B gal = 3 /zG, and a ic / K N ~ 0.45 for 10 /iG. These values are in 
agreement with the expected a ic / K ^ — p — 1 fo r the appropriate ele ctron spectral index p — 1.4 



Moderski et al. 



2005). At e ic > 10 6 , the inverse- 



and 1.6, respectively (see a wide discussion in 
Compton emission breaks to oc e^ 2,6 , again as expected in the case of a steep power law electron 
continuum oc 7~ 3 - 5 at 7 > 7 Dr - In addition, Figure 6 presents the 0.05 — 1000 MeV spectrum of 
the Cen A source as constrained by d ifferent instruments on board CGRO (1991-1995; 'low' and 



'intermediate' states; 



Steinleetal 



1998). As shown, the low-energy (1 — 100 GeV) flat-spectrum 



part of the predicted halo, if detected by the GLAST instrument in the future, could be spectraly 
distinguished from the (presumably nuclear) component observed by CGRO. 

Figure 6 indicates that almost all the energy injected to the elliptical host via annihilation of 
the nuclear 7-ray emission on the starlight photon field is re-emitted via the synchrotron emission 
at ~ 10 13 — 10 14 Hz frequencies, and via the inverse-Compton emission at > 0.1 TeV photon 
energies. Unfortunatelly, the secondary synchrotron photons have almost the same energy as the 
target starlight photons, and at the same time much lower total luminosity, and therefore are not 
likely to be directly observed. However, the secondary 7-ray emission, although also relatively 
weak, is more promising for the detection. In particular, the analysis presented above indicates 
that at high photon energies (> 0.1 TeV) one should expect the photon flux from the galactic pair 
halo 



F- ~ 

± ISO 



4vr (B e 



<io 



-11 



Anj 



ph cm 



(39) 



L cu V 10 42 erg/ S/ 

with the Cen A distance d L = 3.4 Mpc and e = 0.2 TeV. Note, that although the starlight 
radiation absorbs and reprocesses only a tiny (~ 1%) fraction of the nuclear 7-ray emisison, the 
small size of the resulting isotropic galactic pair halos (~ 100 77,) make it quite easy to detect such 
structures. In fact, in the the case of Cen A radio galaxy 100 77, corresponds to roughly 4 kpc, 
or about 4 arcmin, which is still within the active nucleus-centered flux integration region of the 
HESS telescope. Th us, the estimate given in equation 21 can be directly compared with the HESS 



upper limits given in 



Aharonian et al. 



(2005). Such comparison implies that the time-averaged 7- 
ray output of the Cen A nucleus is indeed L in j < 10 42 erg s -1 , which is already a meaningful result. 
This demonstrates an important aspect of the presented analysis: observations of nearby AGNs at 
TeV photon energies can provide important constraints on the time-averaged VHE 7-ray fluxes 
produced in their active centers, even if these sources are not belonging to the blazar sub-class 
(i.e., even if the nuclear portions of their jets are not inclined at small angles to the line of sight, 
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Figure 7. The expected photon fluxes of the Cen A halo for £o = 0.1 GeV (thin lines) and eq = 0.2 TeV (thick lines), corresponding to B = 3 
and 10 fiG (solid and dashed lines, respectively), as functions of the injection luminosity L; n j. Horizontal dotted lines indicates HESS upper limit 
as given in Aharonian et al. 1 2005), and the expected GLAST sensitivity (~ 1.5 X 10 — 9 ph cm -2 s -1 for one-year all-sky survey; Bloom 1996). 



and therefore their direct 7-ray emission is Doppler- hidden). Figure 7 shows in more details the 
expected 0.19 — 5 TeV and 0.1 — 300 GeV photon fluxes of the inverse-Compton emission for the 
discussed Cen A galactic pair halo, together with the indicated HESS upper li mit and the e xpected 



GLAST sensitivity (~ 1.5 x 10 ph cm s for one-year all-sky survey; 



Bloo 



m 



1996). As is 



shown, future observations by HESS and GLAST will be able to put independent constraints on 
both the low- and high-energy spectral parts of the predicted halo. 



4 DISCUSSION AND CONCLUSIONS 

Several BL Lac objects are confirmed sources of variable and strongly Doppler-boosted TeV emis- 
sion produced in the nuclear portions of their relativistic jets. It is more than probable, that also 
many of the FR I radio galaxies, believed to be the parent population of BL Lacs, are TeV sources, 
for which strongly Doppler-hidden nuclear 7-ray radiation may be only too weak to be directly 
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observed (although see 



Aharonian et al. 



2003; 



Beilicke et al. 



2005 



for the case of nearby FR I 



radio galaxy M 87 detected recently at TeV photon energies by HEGRA and HESS Cherenokov 
telescopes). Here we show, however, that about one percent of the total, time-averaged TeV radi- 
ation produced by the active nuclei of low-power FR I radio sources is inevitably absorbed and 
re-processed by the photon-photon annihilation on the starlight photon field, and the following 
synchrotron/inverse-Compton emission of the created and quickly isotropized electron-positron 
pairs. Such a re-processed isotropic radiation could be detected in the cases of at least a few 
nearby FR I radio galaxies, providing interesting constraints on the unknown parameters of the 
active nucleus and the elliptical host. 

In the case of the Cen A radio galaxy considered in this paper, we found that the discussed 
mechanism can give distinctive radiative features due to the isotropic 7-ray emission of the electron- 
positron pairs injected by the absorption process into the interstellar medium of the elliptical host 
(its inner parts in particular, roughly within the radius of 4 kpc from the galactic center), and 
inverse-Compton upscattering thereby starlight radiation to the ^ TeV photon energy range. The 
resulting 7-ray halo is expected to possess a spectral peak at ~ 0.1 TeV photon energies, preceded 
by a flat continuum due to the dominant Klein-Nishina cooling of the radiating electrons, and fol- 
lowed by a steep power-law oc e^J- Tl+0 ' 5 \ where I\ is the photon index of the primary (nuclear) 
7-ray emission. Such a halo should be strong enough to be detected and mapped by stereoscopic 
systems of Cherenkov telescopes like HESS, and, at lower photon energies, by GLAST. 

All of the above findings should apply as well to the other nearby FR I sources. We note in this 
context, that the kinetic power of the Cen A jet, L s ~ 10 43 erg s _1 , is rather low w hen compared 



to ot her FR I sources (e.g., > 10 44 erg s 1 in the case of M 87 radio galaxy: IStawarz et al 



2006b, and references therein). Thus, other (though more distant) objects of the FR I type may 
posses more luminous isotropic halos than Cen A analyzed here. Indeed, taking the time- and angle 
average nuclear 7-ray (~ TeV) luminosity L in j ~ 7] ra _dLj, where r] ra d is the radiative efficiency, 
and the re-processed (> 0.1 TeV) luminosity L iso ~ r L in j with r ~ 0.01 as estimated in this 
paper, one can find that modern Cherenkov telescopes with the available sensitivity limit 10 -13 
erg cm~ 2 s _1 will be able to detect the discussed halos from the objects located within the radius 



rougly d L ^ 100 \/fr? ra ri /0.1) fLj/10 44 er gs -1 ! Mpc. Is it therefo re possible to attribute the 7-ray 



flux detected recently from M 87 system (|Aharonian et al 



2004) to the isotropic halo of its host 



galaxy? The answ er is negative, since v ariability of this emission established on the time-scale of 
months and years ( Beilicke et all2005 ) excludes any extended 7-ray emission sites. 
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Figure Al. Optical depth for photon-photo annihilation due to the dust lane (solid line) and the elliptical host (dotted line). 



APPENDIX A: OPACITY DUE TO THE 'DUST LANE' 

The continuum far-infrar ed emission of NGC 5128, produc ed most likely by massive young stars 



and diffuse cirrus clouds (Jo y et al.ll9 88: 



EckarLsLai 



1990), is concentrated within the 'dust lane' . 



We model this feature as a thin disc, perpen dicular to t he jet axis, centered on the active nucleus, 
and extending up to radii Rr v ~ 4 kpc (see 



Israel 



1998). We also restrict the analysis to A 



fir 



100 



yum radiation, for which we take the total observed (IRAS) flux ~ 400 Jy (iGolombe k et al.1 11988). 
corresponding to the luminosity L fir ~ 1.6 x 10 43 erg s -1 . The number density of the far-infrared 
photons, assumed to be uniformly distributed within the 'dust-lane', is then 

Lsx 



nflr, e (£,fi) 



8 7T 2 RL e fir m e c 3 



5(e- e fir ) 



(Al) 



The appropriate optical depth for the photon-photon annihilation, analogous to the one given by 
equation 13, can be then evaluated as 

-^fir ?"b 



re- 



8 7r 2 Rl. e fir m e c 3 



St 



+i 



dw (1 — zd) <Xw(e 



cos[arccot(g/100)] 



77V<_7> 



(A2) 
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This optical depth is shown on Figure 8 (solid line), together with the optical depth due to the 
starlight emission of the elliptical host evaluated peviously (dotted line). As shown, opacity due to 
the dust lane is much smaller than the opacity due to the elliptical host (because of the different 
spatial distribution of the target photons), and in addition regards only very high energies (> 100 
TeV) of the primary 7-rays. 
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